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Abstract —Integrating power transfer into wireless communica¬ 
tions for supporting simultaneous wireless information and power 
transfer (SWIPT) is a promising technique in energy-constrained 
wireless networks. While most existing work on SWIPT focuses 
on capacity-energy characterization, the benefits of cooperative 
transmission for SWIPT are much less investigated. In this 
paper, we consider SWIPT in an orthogonal frequency-division 
multiplexing (OFDM) relaying system, where a source node 
transfers information and a fraction of power simultaneously 
to a relay node, and the relay node uses the harvested power 
from the source node to forward the source information to the 
destination. To support the simultaneous information and energy 
cooperation, we first propose a transmission protocol assuming 
that the direct link between the source and destination does 
not exist, namely power splitting (PS) relaying protocol, where 
the relay node splits the received signal power in the first hop 
into two separate parts, one for information decoding and the 
other for energy harvesting. Then, we consider the case that the 
direct link between the source and destination is available, and 
the transmission mode adaptation (TMA) protocol is proposed, 
where the transmission can be completed by cooperative mode 
and direct mode simultaneously (over different subcarriers). In 
direct mode, when the source transmits signal to the destination, 
the destination receives the signal as information and the relay 
node concurrently receives the signal for energy harvesting. Joint 
resource allocation problems are formulated to maximize the 
system throughput. By using the Lagrangian dual method, we 
develop efficient algorithms to solve the nonconvex optimization 
problems. 

Index Terms —Simultaneous wireless information and power 
transfer (SWIPT), orthogonal frequency-division multiplexing 
(OFDM), cooperative relay, resource allocation. 

I. Introduction 
A. Motivation and Related Work 

Recently, simultaneous wireless information and power 
transfer (SWIPT) has become an emerging solution to pro¬ 
long the lifetime of energy-constrained wireless networks. 
SWIPT refers to using the same wireless electromagnetic 
wave to transport energy for both information decoding and 
energy harvesting at the receiver. With SWIPT, the integrated 
transceiver designs can be exploited to efficiently use the 
available wireless resources. 

In 0 , 0 , fundamental performance limits of SWIPT 
were studied from an information-theoretic perspective, and 
the tradeoff between the channel capacity and the harvested 
energy at the receiver was characterized. Such capacity-energy 
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tradeoff was also derived for a multiuser system in 0 . These 
papers assume that the received energy can be still harvested 
after passing through an information decoder, which is not 
realizable yet due to practical electronic circuit limitations. 
A practical SWIPT receiver was designed in a, where the 
received signal is split into two separate circuits, one for 
information decoding and the other for energy harvesting. 
The works in a and m studied SWIPT in multiple-input- 
multiple-output (MIMO) channels with perfect and imperfect 
channel state information (CSI), respectively. In particular, the 
authors in 0 proposed two transmission protocols, namely 
“time switching” where the receiver switches between decod¬ 
ing information and harvesting energy, and “power splitting” 
where the receiver splits the signal power into two parts for 
decoding information and harvesting energy. Note that the 
time switching protocol cannot operate information and power 
transfer simultaneously. The time switching protocol and the 
power splitting protocol with power allocation were investi¬ 
gated in 0 and 0 for capacity-energy tradeoff, respectively. 

A few attempts have been made very recently to study 
wireless information and power transfer in cooperative relay 
systems a go). To be specific, the outage probability and 
the ergodic capacity for time switching and power splitting 
protocols were derived in (9). The authors in d considered 
an energy cooperation scenario where each user shares a 
portion of the harvested energy with the other users, and the 
optimization of energy arrivals for throughput maximization 
using Lagrangian duality was developed. Note that lITOll con¬ 
sidered a full-duplex relay which can transmit and receive 
data/energy simultaneously. In 1 11I I. the authors investigated 
wireless information and power transfer in cooperative net¬ 
works where the randomly located relays assist one source- 
destination pair, and outage probability and diversity gain 
were characterized by stochastic geometry. The authors in 
G2l studied the outage performance of power strategies at an 
energy harvesting relay for multiple source-destination pairs. 
The optimal time-switching ratio was investigated in full- 
duplex relaying system under different communication modes 

in Cl- 

On the other hand, orthogonal frequency-division multi¬ 
plexing (OFDM) based relaying transmission is a powerful 
tool to enable high date rates and has been adopted as a 
key technique for the next generation communications. The 
optimization and resource allocation schemes were proposed 
for various settings in OFDM relay systems, e.g., single-user 
single-relay m-m, multi-user single relay CD, and multi¬ 
user multi-relay G8l - f23l . For OFDM-based SWIPT, power 
optimization in OFDMA systems for different configurations 
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was studied in [[24). The authors in 1251 studied SWIPT in 
downlink OFDMA systems with power splitting receivers, 
and suboptimal iterative algorithms were developed to solve 
the non-convex resource allocation problems. In 1261 . the 
weighted sum-rate maximization problem subject to a mini¬ 
mum harvested energy constraint on each user was studied, 
under either time switching or power splitting protocol to 
coordinate energy harvesting and information decoding. 

The above OFDM-based SWIPT Il24l - ll26l focus on tradi¬ 
tional direct transmission. The authors in l27l studied SWIPT 
in MIMO-OFDM relay networks and solved the non-convex 
resource allocation problems in a suboptimal step-wise way. 
It is thus observed that SWIPT in OFDM relay systems is not 
well studied. Using the architecture of OFDM relaying, the 
potential of SWIPT can be fully explored in various dimen¬ 
sions, including frequency diversity and cooperative diversity 
by designing efficient resource allocation schemes. It is thus of 
great importance and necessity to capture the full potential of 
OFDM relaying for SWIPT. This motivates our study. In this 
paper, we aim to answer the following questions: What is the 
transmission protocol to be used for implementing SWIPT in 
an existing OFDM relay network? How to efficiently allocate 
wireless resources, such as power and subcarriers to maximize 
the throughput given the transmission protocols? 

B. Contributions 

In this paper, we consider a three-node relaying scenario, 
where a source node transfers a portion of its energy to a 
decode-and-forward (DF) relay node in return for its assistance 
in information transmission using OFDM. The main contribu¬ 
tions and results of this paper are summarized as follows: 

• Assuming the direct link between the source and des¬ 
tination is not available, power splitting relaying (PS) 
protocol for supporting SWIPT in OFDM relaying is 
considered, where the relay node splits the received signal 
power in the first hop into two separate parts, one for 
information decoding and the other for energy harvesting. 
Then considering the case that the direct link between 
the source and destination is available, we propose the 
transmission mode adaptation (TMA) protocol to sup¬ 
port SWIPT in OFDM relaying, where the end-to-end 
transmission can be completed through relay (coopera¬ 
tive mode) and direct source-to-destination transmission 
(direct mode) simultaneously. For direct mode, when the 
source transmits signal to the destination, the destination 
receives the signal as information and the relay node 
concurrently receives the signal for energy harvesting. 
For both protocols, the relay node harvests power from 
the source node as the transmit power for information 
forwarding, which can be regarded as simultaneous two- 
level cooperation, i.e., information-level cooperation and 
energy-level cooperation, thanks to the parallel structure 
of OFDM relaying. 

• For the PS protocol, the joint problem of power allo¬ 
cation and determining power splitting ratio is studied. 
For the TMA protocol, we consider a joint problem of 
mode selection (between cooperative and direct modes). 
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Fig. 1. (a) Illustration of the PS relaying protocol, (b) Diagram of the relay 

receiver over all subcarriers. 


subcarrier and power allocation. By using the Lagrangian 
duality method, we develop efficient algorithms to solve 
these nonconvex problems. 

« Simulation results show that for the PS protocol, more 
resources are allocated to information transfer and less 
resources for power transfer. The power allocated to 
information transfer is increasing when the relay node is 
closer to the destination; on the other hand, the power 
allocated to power transfer is decreasing in this case. 
For TMA protocol, more power should be allocated to 
cooperative mode when the relay node is located at the 
middle, and less power is located when the relay node 
is closer to either the source or destination. The case is 
reversed for the direct mode. We also show that for both 
protocols, it is better to place the relay node closer to the 
source node rather than in the middle between the source 
and the destination. 


C. Organization 

The remainder of this paper is organized as follows. Section 
II describes the system model and the proposed transmission 
protocols. Section III and Section IV present the joint resource 
allocation problems for the proposed PS and TMA protocols, 
respectively. Section V provides extensive simulations. Finally, 
we conclude the paper in Section VI. 

II. System Model and Transmission Protocols 

We consider a three-node relaying system, consisting of a 
source node S, a relay node R, and a destination node D 
as shown in Figs. fflE OFDM is used at the physical layer 
for the two-hop transmission. The relay node is half-duplex, 
i.e., it cannot receive and transmit simultaneously. The end-to- 
end transmission is completed by two equal time slots. In the 
first time slot, the source node transmits signals over OFDM 
channels while the relay node receives. In the second time 
slot, the relay node forwards the processed information to the 
destination using the DF strategy. 

Using the parallel structure of OFDM relaying, we first 
propose the PS transmission protocol to support SWIPT as¬ 
suming that the direct link between the source and destination 
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Fig. 2. Illustration of the TMA relaying protocol, (a) Cooperative information 
transmission mode, (b) Direct information transmission mode, (c) Diagram of 
the relay receiver. 

does not existQ Fig. □ depicts the proposed PS protocol. 
The signal power on all subcarriers is split into two parts at 
the relay receiver, one for information decoding and one for 
energy harvesting. Specifically, the relay receiver determines 
the optimal power splitting ratio a allocated to the information 
decoder and 1 — a to the energy harvester. That is, all 
subcarriers have the same power splitting ratio. The energy 
harvested during the first time slot is used as a source of 
power for the relay to forward the source information to the 
destination. 

Then, we assume that the direct link between the source 
and destination is available and propose the TMA protocol in 
Fig. [2] By exploring the direct link between the source and 
destination, the end-to-end information transmission can be 
completed by two transmission modes simultaneously (over 
different subcarriers), e.g., cooperative mode (see Fig. |2|a)) 
and direct mode (see Fig. |2jb)). Correspondingly, the OFDM 
subcarriers in the first time slot are partitioned into two groups 
for the two transmission modes respectively. In particular, 
for the direct mode, when the source transmits signal to the 
destination, the destination receives the signal as information 
and the relay node concurrently receives the signal for energy 
harvesting. The harvested energy at the relay node at the 
subcarriers with direct mode is used as a source of power 
to forward the source information for cooperative mode. At 
the relay node (see Fig. [2c)), the signal power from the 
subcarriers with cooperative mode and direct mode are sent 
to information receiver and energy receiver, respectively. 

There are two points about the protocols to be noted. First, 
for both protocols, the harvested powers from all subcarriers 
during the first time slot are added up as a total power 
for information forwarding over all subcarriers in the second 

'This is due to the impairments such as wireless fading, shadowing, path 
loss, and obstacles, which make the effects of direct link negligible. 


time slot. Second, we assume that the information from one 
group of subcarriers in the first hop can be decoded and re¬ 
encoded jointly and then transmitted over a different group of 
subcarriers in the next hop. This is optimal for DF relaying 
over multi-channels |[28l . 

It is noticed that for the PS protocol, we let the relay 
receiver split the signal energy into two parts for all subcarriers 
with the same power splitting ratio, since the power splitting 
is performed in analog domain before the digital OFDM 
demodulation. Moreover, for the TMA protocol, we assume 
that the relay receiver has ideal bandpass filters so that it is 
able to tap into different subcarriers. This is also possible 
in OFDM-based ultra wide band (UWB) system where the 
ratio of bandwidth and carrier frequency is more than 20% 
(according FCC’s definition), which makes the bandpass filters 
not sharp and thus feasible. 

The fading is modeled by large-scale path loss, shadow¬ 
ing, and small-scale frequency-selective Rayleigh fading. The 
transmission from the source to destination is divided into 
consecutive frames. It is assumed that the fading remains 
unchanged within each transmission frame but varies from one 
frame to another. We also assume that perfect CSI is available 
at all nodes. Without loss of generality, the additive noises 
at all nodes are independent circular symmetric complex 
Gaussian random variables, each having zero mean and unit 
variance. 

In this paper, we assume that energy loss and nonlinearity, 
due to processing and circuitry at the transmit/receive, are 
neglected for simplicity. The assumption is widely adopted 
in the prior work PIl- lfToll . li26l . 

III. Power Splitting Protocol 

In this section, we study the PS protocol and the corre¬ 
sponding joint optimization problem of power allocation and 
power splitting. 

A. Problem Formulation 

Let p Sin denote the source power on subcarrier n in the 
first time slot. As stated in l26l . all subcarriers would have the 
same power splitting ratio in an OFDM-based SWIPT system. 
Thus we denote a as the fraction of power for information 
decoding and the remaining 1 — a is the fraction of power 
for energy harvesting. The harvested power from the source 
node over all subcarriers in the first time slot is used as the 
source of transmit power to forward the source information 
to the destination. Let h n and g n are the channel gains over 
subcarrier n in the first and second hops, respectively. Denote 
J\f := {1, ■ ■ ■ , N } as the set of subcarriers, then the transmit 
power of the source {p s ,n} and the relay {p r , n } should satisfy: 

^ Ps,n < P, (1) 

nGAf 

and 

^ ' Pr,n — ^ ' Ps,n^In(l (2) 

new new" 
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In addition, we consider the peak power constraint on each 
subcarrietQ 


PpS ■■= {0 < p s ,n < P n , 0 < p rjU < p n , Wl}. (3) 

To avoid trivial solution, we assume that YneAfP™ — P- 
We can easily obtain the end-to-end achievable rate of the 
PS protocol as 

Rps = ^ min { 22 log(l + p Syn h n a), 22 log(l + p r ,nSVi)}- 

n£Af n£j\f 

(4) 


In the PS protocol, our goal is to maximize the system 
throughput by jointly determining the transmit power and 
power splitting ratio on each subcarrier. Let p = {p s ,n,Pr,n}, 
and the problem can be formulated as 

PI : max Rps (5a) 

{p,a} 

s.t. {j}-©. (5b) 

B. Proposed Solution 

The problem in PI is a mixed integer optimization problem 
and nonconvex. Finding the optimal solution needs exhaustive 
search with prohibitive complexity. However, as shown in 
(29l . for the nonconvex optimization problems in multicarrier 
systems, the duality gap is zero under the “time-sharing” 
condition, and the time-sharing condition is always satisfied 
if the number of subcarriers is large. Therefore, in this paper 
we solve the resource allocation problems in dual domain. 

We first introduce non-negative Lagrangian multipliers //, v, 
A, and k rated with the constraints Q-Q, respectively. The 
Lagrangian of PI is 


£ps(A, K,p, v,p, a) = \ 


1 

2 


+K 



■ 22 l°g(l + Ps,nh n a) 

nGAf 

log(l + Pr,n9n) ~ RpS 


- -Rps 

+ Rps 


+p 

P 'y 1 Ps,n 

+ V 

22 Ps.nM 1 - a ) 

^ ' Pr,n 


nGAT 


.n£Af 

n£Af 


Define V\ as the set of all primal variables that satisfy the 
constraints, then the dual function of PI can be expressed as 


5 ps(A, K,p, v) = max £ P s(A, k, p, v,p, a). (7) 

{p,a}£X>l 

Computing the dual function .gps(A, k, p, v) needs to find 
the optimal {p*,a*} for the given dual variables {A, k, p, v}. 
In the following, we detail the derivations. 

We first investigate the part of the dual function with respect 
to the rate variable Rps- 


5o(A, k) = max (1 - A - k)R p s- (8) 

Rps 


To ensure the dual function to be bounded, we have 1 A — 
k = 0 such that <?o(A, k) = 0, which means that k = 1 — A. It 
is nontrivial that 0 < A < 1 making k non-negative. 

By substituting the result to ©, the Lagrangian becomes 


£ps(A ,p,v,p,a) = -^'22 lo g(l +Ps, n h n a) 


i&M 


1 — A 


H- 7T~ 5 Z lo g(! +Pr,n9n) + M 




P ^ Ps,n 
neAf 


22 Ps,nh n ( 1 - a) - 22 Pr ’ r 
n£Af nEAf 


(9) 


Note that when tackling the problem, the power splitting 
ratio a appears in the achievable rate expression and couples 
the power allocation on each subcarrier, which results in a 
nonconvex problem. To this end, we first derive the optimal 
power allocations for given a in the following proposition. 

Proposition 1: For a given power splitting ratio a, the 
optimal power allocations are 


Pn : 





if vh n { 1 — a) > p 
otherwise. 


( 10 ) 


P 


* 

r,n 


l-A 
2 v 


( 11 ) 


Proof: Note that for a fixed a, the Lagrangian is concave 
in the power allocation variables. For each p s , n , the corre¬ 
sponding sub-Lagrangian is 

Lps(p s ,n ) = ^\og(l+p s , n h n a)-p Stn {p-vh n [l-a)). (12) 

Since A > 0, Lps(Ps.n) is an increasing function of p StU if 
vh n { 1 — a) > p and thus the optimal solution is p* n =p n . 
Otherwise Lps(Ps.n) is a concave function of p s ,n- By the 
Karush-Kuhn-Tucker (KKT) conditions (30), the maximum is 
achieved at £>Z 'g^ Pa,n ' ) = 0. p* n directly follows the KKT 
conditions. ■ 

Substituting the optimal power allocations into the La¬ 
grangian, it is a function of the power splitting ratio a, which 
still turns out to be neither a convex nor a concave function. As 
the closed-form solution of the power splitting ratio a seems 
infeasible, we resort to a numerical method, i.e., we enumerate 
all possible values of a in the interval [0,1] with a sufficiently 
small step-size and select the value of a* that achieves the 
maximum of the Lagrangian. 

After finding the optimal {p*,a*}, we turn to the dual 
problem, which can be expressed as 


min g PS (X 1 p,p) (13a) 

{A.M.i'} 

s.t. 0 < A < l,p > 0,u > 0. (13b) 


2 This is known as “spectral mask” in OFDM systems introduced by FCC, .. , , ,■ ■ ■ 

for considering practical power-amplifier limits and controlling interference Since a dual function IS al ways convex by definition, we 
to other devices. employ the ellipsoid method lf3Ql to simultaneously update 
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Algorithm 1 Proposed Algorithm for PI 

1: initialize {X,p, v} as non-negative values. 

2 : repeat 

3: Compute the optimal power allocations p*(a, A, p,v) 

using Proposition Q] 

4: Find the optimal power splitting ratio ct* (A, p, v) using 

numerical search with step-size e. Here 0 < e < 1 
controls accuracy. 

5: Update {A, p, v) by the ellipsoid method using the 

subgradients defined in (IT4l >- (IT6t . 

6: until (A,/z, i/} converge. 


the dual variables {X 7 p, v} toward the optimal {A*,p* ,v*} 
by using the following subgradients. 

= \ J2 togi 1 + P*s,n h na*) l °g( l +Pr,n9n), 

n£Af n£ J\T 



(14) 

= P - y P*s,n , 

(15) 

n^M 


y P*s,nhn{l - a*) - y p* >n . 

(16) 


n£Af n€LAf 


C. Algorithm and Complexity 

We summarize the proposed solution for PI in Algorithm 
1. The computational complexity of the ellipsoid method is 
0(q 2 ), where q is the number of the dual variables, and 
clearly q = 3 in our paper. Denote e as the step-size of the 
search procedure for the optimal a*, the total complexity of 
the proposed algorithm is 0(q 2 /e). 

IV. Mode Adaptation Protocol 

In this section, we study the TMA protocol and the cor¬ 
responding joint optimization problem of transmission mode 
selection, subcarrier assignment, and power allocation. 

A. Problem Formulation 

Let psr.r, and p s d, n denote the transmit power of the source- 
to-relay and source-to-destination links on subcarrier n at the 
first time slot, respectively. f n represents the channel gain of 
the source-to-destination link on subcarrier n. Binary variables 
y c „ and yd >n are used f° r subcarrier assignment at the first 
time slot, y c>n = 1 indicates that subcarrier n is used for 
cooperative mode and y C}U = 0 otherwise; y d ^ n = 1 indicates 
that subcarrier n is used for direct mode and y,i, n = 0 
otherwise. y c>n and y d , n satisfy the following constraint 

Uc,n + yd,n < 1, Vn € AT. (17) 

Then, the transmit power of the source should satisfy the total 
power constraint: 

(V°.»P sr,n + yd,nPsd,n) — P- (18) 

n£A/" 

The harvested power at the relay node on subcarrier n dur¬ 
ing the direct mode is yd,nPsd,nh n . The total harvested power 


is used as a source of power {p r n } to forward information of 
the source node and satisfy 

^ ' Pr 7 n ^ y ' yd,nPsd,nhn- (19) 

n£A/" n£A/" 

The power allocations are also constrained by 

PTMA •— {0 P Psr,n 5) Pni 9 — Psd,n — Pni 

0 <Pr,n<P n yn}. (20) 

The end-to-end transmission rate of the TMA protocol is 
the sum rates of the direct and cooperative modes: 

Rtma = Rd + R c , (21) 

where 

Rd = y ^ Vd,n log(l T Psd,nfn')i (22) 

n,£A/” 

and 

R c = \ min { 51 Vc ’ n 1 °g( 1 +Psr,nh n ), 

nGAT 

T, y C ,n log(l + Pr,n9n + Psr,nfn) } ■ (23) 

nCAT 

Note that in (l23l i. we assume that the destination can receive 
the information from the source to the relay for maximal-ratio 
combining. 

In the TMA protocol, the goal is to maximize the sys¬ 
tem throughput by jointly allocating power, subcarriers, and 
transmission modes. Denote p = {p S r,n,Psd,mPr,n} and 
V = {yc,m Vd,n}, this problem can be formulated as 

P2 : max 7 ?tma (24a) 

{p,y} 

s.t. O - (ID. (24b) 


B. Proposed Solution 

The problem in P2 is a mixed integer optimization problem 
and nonconvex. To make it more tractable, we first relax the 
binary variables to real ones, i.e., 0 ^ y < 1, and defining new 
variables t = t’ c n , t d , n } with t c>n = y c ,nPsr,n, t' c n = 

yc,nPr,n and t dl n = yd,nPsd,n- 

The continuous relaxation makes t the time sharing factors. 
Moreover, after the relaxation, P2 is a convex problem. We 
notice that applying conventional software packages directly 
on P2 is not sufficient to solve the original problem P2, since 
there is no guarantee that they can return the binary y. 

We can write the Lagrangian as 


£ tma (/ 3 , 7 , < 5 , p , y,t)=y y d}7l log (l + td ’ n ^ n \ 

V Vd,n J 


-X y y c ,n log ( 1 




ga f 


1 — /5 i (-t , t cn 9n + tc,nfn 

lo s 1 


Vc,r 

+ 


■leN 


Vc,n 


P — y ' ( tc,n + td,n) 
n£,V 


+ <5 


'y, td,nh n t c ri 

n&Af n£jy 


( 25 ) 
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where /?, 7, and <5 are non-negative Lagrangian multipliers 
relating to the constraints (OH . ( fl 8 ] >. and (fl9] >. respectively. 
Moreover, the constraint 0 < /3 < 1 must hold. The dual 
function can be expressed as 


3tma(/ 3, T, <5) = max £tma(/3,7, S,p, y, t), (26) 

{p,y,t}eT>2 


where V 2 as the set of all primal variables that satisfy the 
constraints. 

We first get the optimal power allocations in the following 
proposition. 

Proposition 2: The optimal power allocations of the relaxed 
P2 are 


_ * _ 

®c,n — Vc,nP S r,n — Uc,r 


P 


2(7 -6) h n _ 


(27) 


^ c,n yc,nP r ,n Z/c, 

and t*d,n = yd,nP* d ,n with 

P*sd,n = 


(1 ~ P)fn _ P*sr,nfn _ 

2 Sg n g n g n q 


7 — 8h n 


Pn 
J 0 


if 7 > Sh n 
otherwise, 


(28) 


(29) 


Proof: (ITT! ) and (f28l) can be easily obtained by applying 
KKT conditions. The Lagrangian with respect to p s d, n can be 
written as 


^- J d,n(.Psd,n) — lt)g(l T Psd,nfn ) Psd,n[pt (30) 

Since Ld, n is an increasing function of p s d, n if 7 < 5h n , 
thus the maximum of Ld, n is at p* d = p n . If 7 > Sh n , 
Ld^n is a concave function of p s d, n and its maximum is at 

dLd.ni'Psd.'n) _ Q _ 

dps d,n 

We are now ready to find the optimal subcarrier assignment 
between cooperative mode and direct mode at the first time 
slot by the following proposition. 

Proposition 3: For each subcarrier n in the first hop, the 
optimal subcarrier assignment is given by 

f Vc,n = 1 and Vd,n = O’ ^ Jc,n Jd,n q ] ) 

\y*,n = 0 and Vd,u = 1,otherwise, 

where 

Jc,n 7^ lo§(l “1“ Psr,nh'n') ^ cj - P S r,nfn Pr,n9n ) 

-P*«-,n'r-Pr,J, (32) 

Jd,n = l0g(l + P* sd ,nfn) ~ P* sd ,n (7 ~ Sh n ). (33) 

Here p* p* sr n and p* d n are obtained in Proposition [2] 
Proof: To maximize the Lagrangian £ tma , the La¬ 
grangian over each subcarrier should be maximized, which 
can be expressed by 


max 

yc,nJc,n yd,nJd,n 

(34) 

y c ,n,yd,n> 0 


s.t. 

yc,n H - yd,n 5: 1 * 

(35) 


Since the objective function has a bounded objective and its 
maximization over y cn and vy r / Ti is a linear programming 
problem, a global optimum can be found at the extreme points 


Algorithm 2 Proposed Algorithm for P2 

l: initialize {/3, 7 , A} as non-negative values. 

2 : repeat 

3: Compute the optimal power allocations p*(/3,7, 6) us¬ 

ing Proposition [2] 

4: Find the optimal subcarrier assignment y*(/3,'y,S) us¬ 

ing dSp . 

5: Update {/3,7, <5} by the ellipsoid method using the 

subgradients defined in (f37l) -(l39l. 

6 : until {/?, 7, 6} converge. 


of the feasible region ED. Thus we can determine the optimal 
binary solutions of y c _ n and yd.n by exhaustive search over 
{0,1}. In other words, we let the maximizer of „ and Jd. n 
be active and the remaining one be inactive in each subcarrier 
n for the exclusive subcarrier assignment. ■ 

It is worth noting that if J Cj „ = Jd,n, there exist infinite 
maximizers of the objective function, but we are only inter¬ 
ested in the optimal solution in binary form to recover the 
primal exclusive subcarrier assignment constraint (fT71 >. The 
binary recovery for the case of J c ra = J,j n may result in that 
the optimal power allocations obtained in dual domain are not 
feasible for the primal problem. However, the probability of 
J cn = Jd,n is zero, since J c n and J r j 7l are the functions 
of h n and /„ which are the independent continuous random 
variables. Hence, we can uniquely recover the optimal binary 
variables y* for given optimal {/?*, 7 *,6*} with probability 
1 . 

After determining the optimal primal variables {p*,y*,t*}, 
the dual problem can be expressed as 

min 3 tma(/ 3,7’^) (36a) 

{£> 7 . <5} 

s.t. 0 < /3 < 1,7 > 0,J > 0. (36b) 


The optimal dual variables {/3*, 7*, <5* } can be found by 
using the ellipsoid method according to the following subgra¬ 
dients: 


A/3 


2 / , ZZc,n 
neA/” 


y*,r, 


l E y*,n log ( 1 + 

?c*n9n+t* c>n fn\ 

(37) 

y*,n ) ’ 

n£ N ' 


A 1 = P~Y, (C,n + Cj,n)> 

nGAf 

(38) 

A<5 = E ^d,nK 

- E C- 

(39) 

n£j\f 

n£j\f 



C. Algorithm and Complexity 

Finally, we present the proposed solution for P2 in Algo¬ 
rithm 2. The complexity in Proposition [D for all subcarriers is 
0(2N). Combining the complexity of dual updating, the total 
computational complexity is 0(2Nq 2 ). 
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V. Simulation Results 

In this section, we conduct comprehensive simulations to 
evaluate the performance of the proposed protocols. 

We set the distance between the source S and the destination 
D is 2, and the relay node R is located in a line between 
the source S and the destination D. The Stanford University 
Interim (SUI)-6 channel model ll32l is employed for generating 
OFDM channels and the path-loss exponent is set as 3.5. The 
standard deviation of lognormal shadowing is set to be 5.8 dB 
and the small-scale fading is modeled by Rayleigh fading. The 
number of the subcarriers is 64. Moreover, we set p n = 5dB, 
Vro. The step-size of Algorithm 1 for searching a* is set to be 
e = 10~ 3 . In this section, all numerical results are obtained 
by averaging 1000 channel realizations. 

A. Performance of Symmetric Relay 

We first introduce the benchmarks as follows. 

1) Performance Upper Bound of Power Splitting (UB-PS): 
An upper bound of the PS protocol can be obtained by 
assuming that the receive of the relay can dynamically 
determine the power splitting ratio on each subcarrier, 
i.e., a n for all n £ M instead of a. Although the scheme 
may not be implemented in practice, it actually provides 
a performance upper bound for the PS protocol. We give 
the derivations in Appendix lAl 

2) Equal Power Splitting (EPS): For the EPS scheme, the 
power splitting ratio is set to be a = and the pure 
power allocation problem is studied using the algorithm 
for PS protocol. The complexity is 0(q 2 ). 

3) Equal Transmission Modes (ETM): For the ETM 
scheme, the half of subcarriers in the first time slot is 
used for cooperative mode and half for direct mode, and 
the pure power allocation problem is studied using the 
algorithm for TMA protocol. The complexity is 0(q 2 ). 

In this subsection, we consider the case when relay is 
located the middle of the source S and the destination I). From 
Fig. [3] we observe that the TMA protocol outperforms the PS 
protocol as expected. The performance of the PS protocol is 
very close to the upper bound. This suggests that splitting all 
subcarriers with an identical power ratio is very efficient. In 
addition, we observe that as SNR increases, the performance of 
all schemes are bounded at about SNR=25dB, due to the peak 
power constraint on each subcarrier. Moreover, the gains over 
the benchmarks illustrate the benefits of adaptive subcarrier 
assignment and power splitting. 

From Fig. 0] it is found that for the TMA protocol, more 
subcarriers are used for both cooperative and direct transmis¬ 
sion modes when SNR increases. More subcarriers are used 
for cooperative mode when SNR is low, and more subcarriers 
are used for direct mode when SNR is greater than about 20 
dB. 

Fig. 0 shows that the power percentages in the proposed 
protocols. For the PS protocol, one observes that the power 
percentage of information transfer is decreasing with the total 
power constraint of the source node, and the power percentage 
of harvested energy at the relay node is decreasing with the 
source power. For the TMA protocol, the power percentage 



Fig. 3. Sum-rate comparison of the two proposed protocols, where the relay 
is located at the middle between the source and destination. 



Fig. 4. Percentages of used subcarriers for the proposed TMA protocol, 
where the source power is fixed as P = 20 dB. 

of cooperative mode is decreasing with the transmit power, 
and that of direct mode is increasing with the transmit power. 
Moreover, the transmit power used for information transfer 
(cooperative mode) is more than for power transfer (direct 
mode) when SNR is lower than about 20 dB, and the case is 
reversed when SNR is greater than about 20 dB. 

B. Impacts of Relay Location 

In this subsection, we consider the impacts of the relay 
location for system performance. Here, we fix the source 
power as P = 20 dB, and let d as the distance of the source 
node and the relay node. 

In Fig. O we compare the system throughput of the pro¬ 
tocols with different relay locations. It is observed that the 
TMA protocol outperforms the PS protocol over a wide range 
of relay locations. In addition, both protocols perform better 
when the relay node is closer to the source node. Specifically, 
the TMA protocol has the best performance when the relay 
























0.55 



10 20 
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Fig. 5. Power percentages of the two protocols, where the source power is Fig. 7. Power percentages of the protocols, where the source power is fixed 
fixed as P = 20 dB. as P = 20 dB. 



cooperation, two transmission protocols were proposed. The 
corresponding nonconvex resource allocation problems within 
the two protocols were solved with efficient algorithms. 

Appendix A 

Upper Bound of the Power Splitting Protocol 

Define 'D :i as the set of all primal variables that satisfy 
the constraints and a = {a n }, then the dual function can 
be expressed as 


3ps(A, p, v)= max £ PS (A, p, v,p, a), 
{ P ,a}ex> 3 


(40) 


where the Lagrangian is 


£ PS (A, /r, u,p, a) = - log(l +p s ,nh n a n ) 




Fig. 6. Sum-rate comparison of the two protocols, where the source power 
is fixed as P = 20 dB. 

node is located at d = 0.3 and the PS protocol has the best 
performance when the relay node at about d = 0.5. 

From Fig. [7] for the PS protocol, we observe that more 
power is consumed by power transfer when the relay location 
is closer to the source node, and more power is consumed 
by information transfer when the relay location is closer to 
the destination node. For TMA protocol, most power is used 
for direct mode when the relay location is closer to the source 
node, and the power percentage is decreasing with the distance 
of source-to-relay link. Moreover, the direct mode consumes 
minimum power when the relay is located at the middle 
between the source and destination. 


+ ~Y~ X + Pr,n9n) + P 

nEAf 


p - X Ps ’ r 

n£Af 


+V 


^ ' Ps,nhn( 1 Qn) ^ ' Pr,r 
nGjV raSAf 


(41) 


Note that in above we have eliminated the rate variable Rps 
by using the similar method as in © and 0 < A < 1. The 
maximization of £pg(A, /./, v. p. a) can be decoupled into 2N 
subproblems, where N subproblems for the first hop and N 
subproblems for the second hop: 

£ps(A ,[i,i',p,a)= X Lps, n + X Lps, n + (42) 

neAf n£ M 

where 


VI. Conclusion 

In this paper, we studied the wireless information and 
energy cooperation in OFDM relaying systems. The key idea 
is that the source node transfers a fraction of power to the relay 
node in return for its assistance of information transmission. To 
support such simultaneous information cooperation and energy 


~^PS,n — 2 P Ps i n^'n^n') Ps,n vh n {\ Ctn)] ? 


tt _ 
-^PS ,n ~ 


1 - A 


log(l +Pr,ng n ) - Vp r ,, 


(43) 

(44) 


We have the following proposition for deriving the optimal 
power allocation p* and optimal power splitting ratio a*. 
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Proposition 4: The optimal solution to the performance 
upper bound is given by 

1) ct* = 0 and p* sn = p n , if A < 2v and vh n > p. 


2) a* = 1 and p* n = 
2vp n h n < A — 2v, 

3) < = 


_A_1_ 

2/i h n 


J 0 


, if A > 2z/ and 


X—2u 


n - 2up n h n and P*s,n = Pn > if ^ > 2i/, i//i„ > p, and 
2 vp n h n > X — 2v. 


For the relay. 


Pr,n 


1 - A 


1 


(45) 


2^ 5«Jo 

Proof: By taking the derivatives of Lf> g n with respect to 


p S} n and a n , respectively, we have 


dL 


PS.ra 


hnCXri 


dps,r 


A 


2 1 T Ps,nhnOir 


- [m - - a n )], 


5L 


PS,n 


A 


Ps,nhn 


OcX n 2 1 T Ps,nhnQr 

Since 0 < a n < 1, we have 


Ps,nh n V. 


(46) 

(47) 


A Ps,nhn 


2 1 T Ps,nh? 


-Ps,nh n U < 


dL 


PS,n 




. Xp s , n h n , 

< --- Ps,nh n V. 


(48) 


We consider the following cases. 

• Case 1: Xps q hn — p sn h n v < 0. In this case, A < 2v 
and Lps,n is a monotonically decreasing function of a n . 
Thus we conclude that a* = 0. Substituting a* = 0 into 
61, we have 


To maximize Tps,n> the optimal power allocation should 
be 


P 


* 

s,n 


Pn , if Vh n > P 
0 , otherwise. 


(54) 


Combining the above results, Proposition 0] is thus proved. 


Take a closer look at Proposition [4] it is found that if a 
subcarrier is used for pure information transfer, i.e., case 2), 
the optimal power allocation is also water-filling. On the other 
hand, if a subcarrier is used for pure or partial power transfer, 
i.e., cases 1) and 3), the optimal power allocation is a constant 
giving the channel gain is larger than or equal to a threshold 
p/v. 

Then, the dual problem can be expressed as 

min 5 ps(A, p, v) (55a) 

s.t. 0 < A < 1, p > 0, v > 0. (55b) 

Similarly, we use the ellipsoid method to find the optimal 
dual variables {A *,p*,i/*} by using the following subgradi¬ 
ents. 

AX =IY1 lQ g + P*a,n h nOL* n ) fog(l + P*, n 9n), 

n£Af n£Af 


A h = P - P °,n-< 

(56) 

(57) 

n£A f 


P*s,n h n( 1 - O - 

(58) 


n£j\f nGAT 


LpS ,n = Ps , n ( vh n ~ p) ■ (49) 

To maximize Lf, s n , the optimal power allocation should 
be 


The complexity of decomposition is 2TV. Combining the 
complexity of the ellipsoid method, the total complexity of 
the algorithm is 0(2Nq 2 ). 


P 


* 

s,n 


p n , if vh n > P 
0 , otherwise. 


(50) 


• Case 2: x ■ ■^£ h J ln —Ps,nh n v > 0. In this case, A > 2v 

and p S:Tl < 2 uh. V • H ence ^ps n ' s a monotonically 
increasing function of a n . This means that a* = 1. Sub¬ 
stituting a* = 1 into (l43l) and applying KKT conditions, 
the optimal power allocation has the form of water-filling: 


Ps,n 


\_ 

2 p 


1 Pn 


(51) 


- p s , n h n v < 0 and 


A Ps,n.h„ 


Combining the above expression with the conditions that 
A > 2v and p Si „ <p n < we obtain 2 vp n h n < 

A — 2v must hold for this case. 

PaSP V A . P=,nhn 

«_ase j. 2 i + PBtn h„ 

Ps,nh n v > 0. In this case, X > 2v and p Sj „ > , 

which yields to 2vp n h n > A — 2v. By letting (l47l > be 
zero, we have 

* X~2v 

O-n — 


C 2vp sn h n 

Substituting a* into d43b . we have 

L PS,n = Ps , n ( vh n - p) + ^ log 


(52) 


(53) 
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